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Abstract. We consider the reconstruction of a function on a fi- 
nite subset of Z 2 if the line sums in certain directions are pre- 
scribed. The real solutions form a linear manifold, its integer 
solutions a grid. First we provide an explicit expression for the 
projection vector from the origin onto the linear solution mani- 
fold in the case of only row and column sums of a finite subset of 
Z 2 . Next we present a method to estimate the maximal distance 
between two binary solutions. Subsequently we deduce an upper 
bound for the distance from any given real solution to the nearest 
integer solution. This enables us to estimate the stability of solu- 
tions. Finally we generalize the first mentioned result to the torus 
case and to the continuous case. 



1. Introduction 

The basic problem of discrete tomography is to reconstruct a function 
/ : A — > B where A is a finite subset of ll and B a finite subset of R, 
when the sums of the function values along the lines in a finite number 
of directions are given. In this paper we consider I = 2. There is a vast 
literature on the special case A = G Z 2 |0 < i < m, < j < n}, 

B = {0, 1} where the problem is to find the function values from the 
given row and column sums. In 1957 Ryser [23] and Gale [17] indepen- 
dently derived necessary and sufficient conditions for the existence of 
a solution in this special case. Ryser also provided a polynomial time 
algorithm for finding such a solution. However, the problem is usu- 
ally highly underdetermined and a large number of solutions may exist 
[25] . Therefore the quest is often to find a solution of a certain type. 
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For some classes of highly-structured images, such as hv-convex poly- 
nominoes (the l's in each row and column are contiguous) polynomial 
time reconstruction algorithms have been developed (see e.g. [5], [TT] . 
[12]). Woeginger [26] presented an overview of classes of polyminoes 
for which it is proved that they can be reconstructed in polynomial 
time or it is proved that reconstruction is NP-hard. Batenburg [6] de- 
veloped an evolutionary algorithm for finding the reconstruction which 
maximises an evolution function and showed that the algorithm can be 
successfully applied to a wide range of evolution functions. 

We consider solutions as vectors with the values of / as entries. 
Hajdu and Tijdeman [TS] observed that the set of binary solutions is 
precisely the set of shortest vector solutions in the set of functions 
/ : A — y TL with the given line sums, provided that such solutions exist. 
They also showed that the solutions / : A — > Z with the given line sums 
form a multidimensional grid on the linear manifold which consists of 
all the solutions / : A — > K. with the given line sums. Moreover, they 
determined the dimension of this manifold and indicated how to find a 
set of generators of the grid. Later they used their analysis to develop 
an algorithm to actually construct solutions / : A — > {0, 1} in [2"U] . 
whereafter Batenburg [7], [S] constructed much faster algorithms. 

For functions / : {0, 1, . . . , m — 1} X {0, 1, . . . , n — 1} — > K and given 
row and column sums we deduce a new and explicit expression for 
the projection vector f from the origin onto the real linear solution 
manifold in Section 3. We do not know a similar expression for other 
sets of line sums, but show the result can be extended to the torus case 
(A = (Z/nZ) 2 ) and to the continuous case (A = [0,m] x [0,n] C M 2 ) 
in Sections 6 and 7, respectively. In these sections we answer questions 
posed by Joost Batenburg. 

In many applications it suffices to find a solution or almost-solution 
which is guaranteed to be similar to the original. If all the solutions are 
similar, then they will also be similar to the original and we say that the 
solution set is stable. Alpers, Gritzmann and Thorens [I] showed that 
already a small change in the data can lead to a dramatic change in the 
image. Research on uniqueness and stability of solutions was carried 
out by Alpers et al. PQ, [3], [2] and Van Dalen [13] , [H] , [15] , [16] in case 
only row sums and column sums are given. Their estimates depend on 
only few parameters. More general situations were studied by Brunetti 
and Daurat [10] and by Gritzmann, Langfeld and Wiegelmann [18] . In 
Section 4 we use the distance estimates for the solutions to derive new 
stability results. They involve more parameters, but, at least in the 
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given example, yield better results than the estimates obtained by Van 
Dalen. 

In Stolk [23] a system of line sums is called compatible if and only if 
a real solution exists. (Then the projection vector / is an example of 
such a real solution.) Stolk showed that if the line sums are integers and 
they are compatible, then there exists an integer solution. In Section 
5 we derive an upper bound for the Euclidean distance from / to the 
nearest integer solution. 

The interest in discrete tomography arose from the study of atom 
positions in a crystal, but the developed theory has also applications 
in medical imaging and in nuclear science, see [211 122] • The results in 
the present paper are merely of theoretical interest, but applications 
of the method can be found in [5] and they may help to estimate how 
many directions are needed to be sure that the solution is unique so 
that one is certain to have found the original configuration. 



We use the same settings as in \±9\. Let a,6 G Z with gcd(a, b) = 1 
and a > 0. We call (a, b) a direction. Put 



By lines with direction (a, b) we mean lines of the form ay = bx + t (t G 
Z) in the (x, y)-plane. 

Let m and n be positive integers and 



Let R be an integral domain such that R[x, y] is a unique factorization 
domain. (We apply the results from this section both with R — R and 
with R — Z.) If g : A — > R is a function, then the line sum of g along 
the line ay = bx + 1 is defined as ^2g(i,j) where the sum is extended 
over all G A with aj — bi + t. 

We often consider a function v defined on A as a vector with the 
mn values of v as entries. If we want to emphasize this, we write v 
instead of v. We always assume that the coordinates of these vectors 
are arranged according to the elements of A in lexicographical order. 

If R C R, the length of v (or v) equals \v\ = \v\ = \/Yl(i j)£A v (^ j) 2 - 



2. Notation and general results 




A := G Z 2 : < i < m, < j < n}. 
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A set of directions S = {(a d , b d )} k d=l is called valid for A if Yld=i a d < 
m and Yld=i 1^1 < n - Suppose S = {(a d , b d )} d=1 is a valid set of direc- 
tions for A. Write M = ^J =1 a d , N = J2 d=1 \b d \. Put 

k 

Fs(x,y) = Ylf(a d ,b d )(x,y) and F^ iVtS) (x, y) = x u y v F s (x,y) 

d=l 

for < u < m — M, < v < n — N . For these values of u and v define 
the functions mi U)V ,s) '■ A — > R by 

™>(!t t v,S)(hj) = coeff. of OV) in for (z,j) e A. 

The mr u ,v,s) ,s are called the switching elements corresponding to the 
direction set S. By the bottom-left corner of the switching element 
m (u,v,s) we mean the lexicographically first element of A for which the 
function value of mr Ul v,s) is nonzero. It is obvious that F( u>Vt s) £ %[x, y)- 
As noticed on page 122 of [TS], it follows from the above definitions that 
the function value of Fr UtVj s) at its bottom-left corner is ±1. 

We refer to the situation as described above as the general case, and 
to the case that k = 2 and the directions are (1, 0), (0, 1) as the simple 
case. In the latter case we only have row and column sums. 

The following result characterizes the structure of the solution set 
over the reals. 

Lemma 2.1 ([H], Theorem 1). Let S be a valid set of directions for 
A. Then any function g : A — >■ R with zero line sums along the lines 
corresponding to S can be uniquely written in the form 

m—l-M n—l-N 

g = /2 ^2 Cuvm (u,v,s) 

u=0 v=0 

with c uv G R. On the other hand, every such function has zero line 
sums along the lines corresponding to S. 

Let A, S and / : A — > R be given. Then we denote by /o the vector 
which is orthogonal to the (m — M)(n — iV)-dimensional linear subspace 
of the vectors g such that g has zero line sums with respect to S and 
for which /o has the same line sums as / with respect to S. 
We shall also use the following consequence of Lemma 2.1. 

Lemma 2.2 (|l9j. Corollary 1). Let K be the set of the bottom-left 
corners of the switching elements of A with respect to S. Then for 
any f : A — >■ R and for any prescribed values for the elements of K , 
there exists a unique g : A — > Z having the prescribed integer values at 
the elements of K and having the same line sums as f along the lines 
corresponding to S. 
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It further follows from Lemma 2.1 that there are exactly MN~Yld=i a d\^d 
linearly independent homogeneous linear dependencies among the line 
sums. 

Let A, S and the line sums along the lines with respect to S be given. 
Checking whether there is a real solution is simple. It suffices to solve 
the corresponding set of linear equations and to see whether it admits 
a solution. If the line sums are integers and the system is solvable, then 
obviously /o is a rational solution. However, it is not obvious that there 
is an integer solution. The existence of such solutions is guaranteed by 
a result of Stolk (cf. p. 20 of [21]). 

Lemma 2.3 (Stolk, |24|, Corollary 3.2.10). Let A, S be given. Suppose 
all the line sums into the directions of S are integers and there exists 
a real function g : A — > K satisfying the line sums. Then there exists a 
function f : A — » Z satisfying the line sums. 

For the convenience of the reader we sketch the proof. 

Proof sketch. Denote the lines in the directions of S by L and the num- 
ber of directions of S by d. The proof proceeds by omitting elements 
from A which are known to have integer values whence it follows that 
the remaining elements have integer sums along the lines of L. 

If a line of L contains only one element, then this element has an 
integer value. Therefore we can omit it from A and the remaining 
elements have integer sums along the remaining lines of S. Denote the 
remaining elements by A and the lines of S containing at least one 
element of this A by L. If the number of remaining directions is less 
than d, we are finished. 

Otherwise every line of L contains at least two elements of A. Con- 
sider the convex hull of A. This polygon contains at least 2d line 
segments, two into each direction of S. By construction the convex 
hull of the switching element also contains two line segments into each 
direction of 5*. Moreover, by the gcd-conditions it is the smallest non- 
trivial polygon with this property. Hence the convex hull of A contains 
a switching element. Consider the bottom- left corners of the switching 
elements contained in the convex hull of A and take the lexicograph- 
ically smallest among them. According to Lemma 2.2 we can choose 
any integer value for this element. We omit this element from A. The 
remaining elements have integer line sums along 5". 

From now on we proceed by iteration until A becomes the empty 
set. If there is a line with only one element of A we omit this element. 
Otherwise we select the first bottom-left corner of a switching element 
contained in A, choose an integer value for it and omit it. 
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3. Explicit expression for the shortest real solution 

in the simple case 

In case of only row and column sums we give an explicit form of / . 
For this we simplify our notation. Let q (i = 0, ... ,m — 1) and rj 
(j — 0, . . . , n — 1) denote the column sums and row sums, respectively. 

n— 1 m—1 

Further, write D — r j- Note that we have D — °i- 

j=0 i=0 

Theorem 3.1. For any G A we have 

Mi,j) = ^ + ^-—. 

n m mn 

Proof. To prove the statement, we need to check two properties: 
fo is a solution, and / is orthogonal to H. 

We start with the first property. Obviously, for any i — 0, . . . , m — 1, 
we have 

n— 1 n— 1 / „ \ j — » _. n—l 



n-± n-± / n n. n 1 " _1 

e/°(m)=£^+^--)=c--+-e 

\ n m mn I mm 

3=0 3=0 V 7 3=0 

Similarly, for any j — 0, . . . , n — 1, 

m— 1 m— 1 / „ \ n 1 "i—l 

£A(u) = E(- + Ii --)='i-- + i E« = 

' z — ' \n m mn I n n z — ' 

i=0 i=0 v 7 i=0 

which confirms the first property. 

To prove the second property we check orthogonality for arbitrary 

— * — * 

h G H. Since for any h G H with the corresponding h : A — > K both 

n—l 

^^,j)=0 (t = 0,...,m-l) 

3=0 



and 



m—1 



EMm)=0 (j = 0,...,n-l), 

i=0 

we have 

m— 1 n—l m—1 n—l , n \ 

£ E /o(m>(^-) = E E ( - + - - — ) h dJ) = 

1=0 3=0 i=0 3=0 v 7 

m—1 / n \ n— 1 n—l m—1 

= £(-- — )£ + E - E Mm-) = o, 

«=0 v 7 j=0 j=0 1=0 

and the theorem follows. □ 
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Remark. For the proof of Theorem 3.1 it seems to be essential that 
every two lines into the same direction contain the same number of 
pixels, i.e. elements of A (cf. Sections 6 and 7). 

4. Binary solutions 

In this section we consider the general binary case / : A — >■ {0, 1}. 
Let D := Yli{ab)&A /( a ' fy- Then D equals the sum of all the line sums 
in an arbitrary direction (ad, bd). We present a method to give a lower 
bound for the number of correct pixel values in an approximate solu- 
tion and to give an upper bound for the number of entries where two 
solutions can be different. Such results have been obtained in the sim- 
ple case by Van Dalen [15] by a completely different method. She calls 
a function F : A — > {0, 1} uniquely determined if there is no other 
function A — > {0, 1} having the same row and column sums. Let a(F) 
be half of the sum of the absolute differences between the row sums 
of F and the row sums of a uniquely determined function Fq with the 
same column sums as F . Then Van Dalen derived upper bounds for 
the number of places where two solutions F\ and F 2 can differ in terms 
of m,n,D = D(F\) and a(Fi). We shall derive some other estimates 
and then compare these results with those in [15] . 

The following result, which is used in the proof of Theorem 2 of [19], 
shows that the Euclidean distance between f and any binary solution 
is fixed. 

Theorem 4.1. For any solution g : A — >■ {0, 1} we have 

\g-fo\ = \/D-\f \>. 

Proof. Observe that if g is a binary solution then \g\ = \/T). This 
means that such solutions are situated on a hypersphere with the ori- 
gin as center, and of radius \J~D. According to Lemma 2.1 the solutions 
g are located on a linear manifold of dimension (m — M)(n — N) orthog- 
onal to Jq. The intersection of this manifold and the hypersphere is 
a hypersphere (of the appropriate dimension) having fo as center. By 
the theorem of Pythagoras we get that the radius of this hypersphere 

is \J D — |/o| 2 ) an d the theorem follows. □ 

Put (x) = min(|x|, |1 — x\) and E = Y^ILi Sj=i(./o(^ j)) 2 - Then the 
Euclidean distance between f and the nearest integer vector with en- 
tries in {0, 1} is exactly \/E~. Hence we have the following consequence 
of Theorem 14.11 
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Corollary 4.1. If E + |/o| 2 > D, then there is no binary solution. 
If E + |/o| 2 — D, then the only solutions g : A — > {0, 1} are obtained 
by rounding each fo{i,j) to the integer or 1 which is nearest to it. 

Note that the only entries where the rounding is not unique are those 
with value 1/2. If such entries do not exist, the solution is unique. 

If D — E — |/o| 2 is positive, but not too large, we still may conclude 
that a certain fraction of the rounded values agrees with any solution 
g : A {0, 1}. (In most cases we cannot tell which rounded values 
are correct and it may even be impossible to do so.) Suppose the 
rounded value F(i,j) G {0,1} of fo(i,j) is not the right value. If 
x := foihj) — 1/2, then, when replacing the value 1 by 0, in E we 
have to replace (foihj)) 2 — (1 — x) 2 by x 2 . Hence the contribution 
increases by 2x — 1 = 2f (i,j) — 1. Similarly, if fo(i,j) < 1/2, then the 
contribution to E increases by 1 — 2x = 1 — 2f (i,j). 

Order the values \2f (i, — in nondecreasing order, 6 1; b 2 , . . . , b mn , 
say. According to Theorem 14.11 D — |/ | 2 equals E plus the sum of the 
values bi which correspond to wrong values in F. Let s be the value 
with 

(1) h + • • • + b s < D - E - |io| 2 <&! + ... + b s+1 . 

Then at most s pixels can have wrong values. Therefore at least mn — s 
pixels have the right values. Similarly, let t be the value with 

(2) &! + ■ ■ • + b t < 2(D - E - \f \ 2 ) <&i + ..- + b t+1 . 

Then for any two solutions at most t corresponding pairs of pixels can 
have different values. Therefore at least mn — t such pixels have the 
same values. So we have derived the following result. 

Theorem 4.2. (a) Let s be defined as above. For any solution g : 
A — > {0,1} we have g(i,j) = F(i,j) for at least mn — s pairs 
(i = l,...,m;j = l,...,n). 

(b) Let t be defined as above. For any two solutions gi,g 2 '■ A — > 
{0,1} we have gi(i,j) = gi^hi) for at least mn — t pairs (i = 

1, ...,m; j = 1, . . .,n). 

This result can be slightly improved if there are line sums over the 
F(i,j) which do not agree with the corresponding sum over the f(i,j). 
Consider some direction with such line sums. We know that some val- 
ues of F have to be wrong and can therefore increase the value of E by 
securing that among bx,...,b s and bi,...,b t in (pQ) and <$fy, respectively, 
there are not too many representatives from some row. This may yield 
smaller values of s and t, hence better results. 
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Example 4.1. Let m = 6,n = 5. Let the row sums be given by 
5, 4, 3, 2, 1 and the column sums by 4, 4, 3, 2, 1, 1, respectively Then 
D = 15 and, according to Theorem 13.11 30/o is given by 

34 34 28 22 16 16 

29 29 23 17 11 11 

24 24 18 12 6 6 

19 19 13 7 1 1 

14 14 8 2 -4 -4 

Thus |/o| 2 — 166/15. Further we get the following table F by rounding 
the elements of fo to the nearest integer. 

111111 
11110 
1110 
1 1 


A calculation gives E = 13/5, hence D — \ fo\ 2 — E = 4/3. This yields 
s = 9,t = 13. Thus, by Theorem I4.2[ every solution differs at most 
at 9 places from F and any two solutions differ from each other at 
most at 13 places. In fact, the solutions below, due to Van Dalen [15], 
show that the actual numbers s and t can be as large as 8 and 10, 
respectively: 

11110 1 111110 
111010 111100 
110100 111000 
101000 110000 
010000 000001 



Theorem 2 of [15J implies that two solutions can differ at most at 
20 places. This estimate is worse than the 13 places obtained above. 
Van Dalen's estimate depends on few parameters, but holds only for 
the simple case. Theorem 4.2 can be applied for any set of directions. 

Since in the above example F does not satisfy the required row sums, 
an improvement of the value of E is possible. The column sums are 
correct, but the top row sum 6 is 1 too high and the bottom row sum 
is 1 too low. For correct line sums only one of both 16 's on the top row 
may be decreased and only one of both 14's on the bottom row may 
be increased. It follows from the corresponding formula (JT|) that every 
solution differs at most at 8 places from F. The above examples show 
that 8 is the best possible. The same reasoning applied to (E} leads to 
an upper bound 11 for t. 
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The example can be generalized as follows. Let n, q be positive 
integers and set m = (n + l)q. Let the row sums be given by Cj = 
n — 1 for 1 < j < q, Ci q+ j = n — I for 1 < I < n, 1 < j < q, and 
Cj = 1 for nq + 1 < j < [n + l)q. In [15] Van Dalen showed that 
t < 2q\J ' Anqin + 1) + 1 — 2q and that t = 2nq can be reached. After 
an elaborate computation the argument in the present paper yields the 
estimates s < 2nq and t < Anq — 2q. Therefore we have obtained 
an improvement by roughly a factor ^fq compared to the estimate for 
t from |15j . Notice that the present bound for t cannot be further 
improved by a factor 2. 

5. The function values are in Z 

In this section we consider the case of line sums for functions 
/ : A — > Z. We derive an upper bound for the solution nearest to 
some given function h : A — > R in the Euclidean sense. The function 
h can be considered as a prescribed model for the integer solution / 
satisfying the line sums. By applying the result with h = f we derive 
an upper bound for the shortest integer solution. 

Let, as before, S = {(ad, bd)}d =1 be a set of directions. Put again 

k M-lN-l 

F s (x,y) = Y[f(a d ,b d ){x,y) = J2J2 Fi ^ iyj - 

d=l i=0 j=0 

Let R(S) = J2tLo 1 Ef=o F tr Obviously R{S) < 2 k , but the following 

examples show that it can be much smaller. 

If S = {(1, 0), (0, l)},then R(S) = 4, 

if S = {(1, 0), (0, 1), (1, 1), (1, -1)}, then R(S) = 8, 

if S = {(1, 0), (0, 1), (1, 1), (1, -1), (2, 1), (2, -1), (1, 2), (1, -2)}, then 

R(S) = 24. 

Theorem 5.1. Let h : A — > R have integer line sums with respect to 
A. Then there exists a function f : A — >• Z with the same line sums 
with respect to A satisfying 

\f- h\ < l -^R{S){m-M){n-N). 

In the proof we shall use the following lemma. 

Lemma 5.1. Let be given a d- dimensional parallelepiped V whose edges 
have lengths h,h, ■ ■ ■ Then the distance from any point of the par- 
allelepiped to the nearest vertex is at most + l\ + • ■ ■ + l 2 d /2. 

Note that the bound is the best possible and that it is attained by the 
centre point of a hyperblock with edge lengths li, l 2 , ■ ■ ■ , Id- 
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Note further that each face Q of V is generated by d — 1 out of the d 
edge directions, and that all the vertices of Q are vertices of V . 

Proof. By induction on d. Let P be a point of the parallelepiped V. 
Let Q be the point on the boundary of V which is nearest to P. Let i 
be the direction which does not occur among the generators of a face 
Q to which Q belongs. Since on the line through P into the direction i 
the point P is in between two boundary points of V at distance U, the 
distance between P and Q is at most U/2. 

We claim that Q is the orthogonal projection of P onto the 
(r— l)-dimensional hyperplane containing Q. Let Q' be the orthogonal 
projection of P onto Q. If Q ^ Q' , then the distance between P and 
Q' is smaller than the distance between P and Q. As Q' is not in V, 
there is a point on the boundary of V which is in between P and Q' . 
This is impossible. Thus Q is the orthogonal projection of P onto the 
face Q of V and the distance between P and Q is at most h/2. By the 
induction hypothesis applied to Q, the distance from Q to the nearest 

vertex of Q is at most ^Jt\ + V l'f_ 1 + / 4 2 +1 + h /^/2. Hence, by 

Pythagoras' theorem, the distance from P to the nearest vertex of Q 
is at most + • • • + l\j2. Since every vertex of Q is a vertex of V, 
the lemma follows. □ 

Proof of Theorem \5.1[ By Lemma 2.3 there exists a function g : A — > Z 
satisfying the same line sums as ft with respect to 5*. It follows from 
Lemma 2.1 that the (m — M)(n — N) vectors in the ran- dimensional 
linear space corresponding to the m by n blocks gi u ,v,S) generate a 
lattice L in the (m — M) (n — iV)-dimensional subspace of vectors which 
correspond to the m by n blocks of integers with all the line sums 
equal to 0. Thus the set of integer solutions is given by g + L. This 
grid generates parallelepipeds with edges of length a/ R(S) and integer 
solution vectors as vertices which cover the entire real solution space. 
Therefore the real solution vector h is in one of these parallelepipeds. It 
follows from Lemma [5. II that the distance from h to its nearest lattice 
point is at most R(S)(m — M)(n — N). Since the nearest lattice 

point corresponds to an integer solution /, this proves the theorem. □ 

If there exists a solution / with all entries in {0, 1}, then this solution 
is the shortest among all integer solutions, since 

m n m n 

^=EE/(i.j)<EE/(u') 8 =i/? 
i=i j=i i=i j=i 

and equality holds if and only ii f : A —¥ {0,1}. If there does not 
exist such a solution, one may ask for the shortest integer solution. 
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The following upper bound for the Euclidean length of the shortest 
integer solution is an improvement of Theorem 2 of [TH] where the 
corresponding upper bound is (2 fe_1 + l)y/mn. 

Corollary 5.1. Let f : A — > M and a set of directions S be given such 
that the line sums of f with respect to S are integers. Let fo be the 
projection vector with respect to S and the line sums. Let g : A — > Z 
be the shortest integer solution in the Euclidean sense having the same 
row and column sums as g has. Then 

\g- fo\ < ^R{S){m-M){n-N). 

Proof. Apply Theorem 15. II with h = fo- This yields a solution 
g : A -»■ Z with \g - f \ < \/R(S) (m -M)(n-N). Since the short- 
est integer solution is the integer solution which is nearest to fo, the 
inequality holds for the shortest solution too. □ 

Lower bounds for \g — / | may be obtained by adapting the methods 
from Section 4 and applying the inequality of Cauchy-Schwarz. 

6. The torus case 

Let n be a positive integer and 

A := G Z 2 : < i < n, 0<j<n}. 

In what follows, by (u,v) G Z 2 we shall mean the unique point {i,j) 
of A with i = u (mod n) and j = v (mod n). (That is, we identify A 
with the torus Z 2 (mod n).) 

A direction (a, b) e Z 2 is admissible for A if < a < n, < b < n 
and gcd(a, b) = 1. Let (a, b) be an admissible direction for A. Then 
the torus lines into the direction (a, b) in A are defined by 

l(a,b),(i,j) = {(x,y) G A : bx-ay = bi-aj (mod n)} G A). 

The next statement shows that torus lines can be obtained from 
ordinary lines of Z 2 . 

Lemma 6.1. Let (x,y) G A, and lt a ,b),(i,j) be a torus line in A. Then 
{x,y) G l(a,b),(i,j) if and only if there exist integers io,jo su °h that x = i 
(mod n), y = j (mod n) and bi — aj = bi — aj . 

Proof. Since x = io (mod n), y = j (mod n) and bi — aj = bi — aj 
imply bx — ay = bi — aj (mod n), the 'if part is obvious. 

To prove the 'only if part, let (x, y) G l( a ,b),(i,j), i-e. bx — ay = t 
(mod n) with t := bi — aj . Then there exists an integer c such that 



(3) 



bx — ay = cn + t. 
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Since gcd(a, b) = 1, there are integers xq and yo such that c = ayo — bxo- 
Substituting this into ([3]), we obtain b(x +nx ) — a(y +ny ) = t. Writing 
i = x + nx and jo = y + ny , the statement follows. □ 

Lemma 6.2. Let lf a ,b),(i,j) be a torus line in A. Then we have 

l(a,b),(iJ) = {(*> J) + t(a,b) : t = 0, 1, . . . ,n - 1}. 

Proof. Obviously, +t(a, b) : t = 0, 1, . . . , n - 1} C /( ,6),(i,i)- Let 

now (a;,y) G A be any point of l< a ,b),{i,j)- Then by Lemma EH there are 
integers i , j with i = x (mod n) and j = V (mod n) such that we 
have 

(4) bi - aj = bi- aj. 

By gcd(a, b) = 1 this implies a \ i — i and b \ jo~j, whence i = i+t±a 
and jo — 3 + for some integers t 1? t 2 . If a6 = then we clearly 
may take t\ = t%. Otherwise, from (j3j) we get 6tia — at?b = implying 
ti =t2- Thus in any case we have 

(x,y) = (i ,jo) = (ij) +t(a,b) (mod n) 
for some t 6 {0, 1, . . . , n — 1}. Hence the statement follows. □ 

Lemma 6.3. Every torus line in A goes through precisely n points of 
A. Further, two distinct torus lines in A into the same direction have 
no points in common. 

Proof. Let l( a ,V),(i,j) be a torus line in A. Noting that for < t\ < t 2 < n, 
the points +h(a, b) and +£2(0, b) are obviously distinct, the 
first part of the statement immediately follows from Lemma 16.21 

To prove the second statement, assume that + ti(a,b) = 

(22, J2) +t 2 (a,b) with some (n, ji), (z 2 , 32) € A and < ti < t 2 < n. 
Then we have = {ii-,3i) + (^2 — ti)(a,b), whence obviously, 

l(a,b),(h,ji) Q ^(a,6),(i 2 ,j2)) which by symmetry implies the second state- 
ment. □ 

In view of Lemma [6T3l for any admissible direction (a, b) for A there 
exist n points Pi (a, b), ... , P n (a, b) in A such that the torus lines with 
direction (a, b) through these points are parallel, and moreover, each 
point of A belongs precisely to one of these lines. 

If (7 : A — > R is a function, then the torus line sums of g along an 
admissible direction (a, b) for A are defined by 

L g ((a,b),t) := 9(hj) (t=l,...,n). 
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Let S be an arbitrary set of admissible directions for A and g be any 
fixed function mapping A into ffiL Then the functions / : A — > R with 

L f ((a, b),t) = L g {{a, b),t) ((a, b) G S, t = 1, . . . , n) 

2 

clearly can be identified with a linear manifold in W 1 . Write H g for 
the set of such functions /. Further, write 

n 

t=i 

for the total sum of g (sum of the line sums in a direction (a, b)). 
Observe that D g is independent of the choice of (a, b). 

Lemma 6.4. Let (u,v) G A and (a,b), (c,d) be admissible direc- 
tions for A with gcd(ad — bc,n) = 1. Then the torus lines l( a ,b),(i,j) an d 
l(c,d),(u,v) have precisely one point in common. 

Proof. Clearly, the common points of the torus lines l( a ,b),(i,j) an d l( c ,d),{u,v) 
belong to the solutions (t, r) of the system of congruences 

ta — rc = u — i (mod n) 
tb — rd = v — j (mod n) 

Since by gcd(ad — be, n) = 1 the above system has precisely one solution 
(t, r) modulo n, the statement follows. □ 

Two admissible directions (a, b) and (c, d) for A are called indepen- 
dent, if gcd(ad — bc,n) = 1. 

Theorem 6.1. Let S = {(a\,bi), . . . , (ak,bk)} be a set of pairwise in- 
dependent admissible directions for A, and let g : A — > R. Then the 
shortest element f G H g (i.e. the one with smallest Euclidean norm) 
is given by 

z — / n n z 

(a,b)eS 

where t(i,j) is that index from {1, . . . ,n} /or which belongs to the 
Une Z( 0l 6),p t(w) ■ 

Proof. To prove the statement, we need to check two properties: that 
/o is a solution, and that / is orthogonal to H , i.e. to the subspace 
H g with g being identically zero. 

We start with the first property. Take any torus line 

l(a,b),p t (a,b) = Pt(a, b) + r(a, b) (t G {1, . . . , n}, r = 0, 1, . . . , n - 1). 
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Further, let (c,d) G S with (c,d) ^ (a,b). Then by Lemma 16.41 each 
torus line P s (c,d) + r(c,d) has precisely one point in common with 
l(a,b),p t (a,b)- By the definition of f this yields, for (a, b) G S, 

J2 fo(P t (a, b) + r(a, b)) = L 9 ((a, b),t) + tlS - (£zi2k 
z — ' ran 

r=l 

which shows that /o G 

To prove the second property, let h G Ho be arbitrary. Then we 
clearly have 

n-1 

^ /i(P t (a, b) + r(a, 6)) = ((a, 6) G S*, t = 1, . . . , n). 

r=0 

Thus 



fo{h3)Kh3)= }^ ~ ^ — W> 

(iJ)eA (ij)eA \(a,6)eS ' / 

= £ f) ( L »» a - ') - gL^S) g t) + r(a, 6)) = 

and the theorem follows. □ 

7. The continuous version in the simple case 

Let m and n be positive real numbers and let A be a Lebesgue- 
measurable subset of T = [0, m] x [0, n]. Write (x, y) for the charac- 
teristic function of A inside T. 

For any bounded Lebesgue-measurable function / : T — > R, let 
Cf(x) : [0, m] — > R and rf(y) : [0, n] — > R be the column integrals and 
row integrals of f(x,y), that is 



o 



c f( x ) = I f( x >y)dy ( x e [o,m]) 

and 

r /(j/) = / f( x ,y) dx (ye[o,n])), 



o 

respectively, where integration is always meant in the Lebesgue sense. 
Since f(x, y) is bounded, by the theorem of Fubini we know that these 
functions exist. Note that the same is true for Cf A (x) and Tf A (y). 

Let L denote the set of bounded Lebesgue-integrable functions T — > 
R having column integrals Cf A (x) (x G [0, m\) and row integrals rf A (y) 
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(y G [0, n\). Further, write H for the set of bounded Lebesgue-integrable 
functions T — > R having vanishing row integrals and column integrals. 
Observe that H is a closed linear subspace of the linear space C of 
bounded integrable functions T — > R. Further, for any gi,g% G L we 
obviously have <?i — #2 G i/. In other words, L = g + H with any g E L. 
We recall that the inner product in £ is given by 

(f(x,y),g(x,y)) = / f(x,y)g(x,y)dxdy 



for /, g G £. The following theorem describes the shortest element in 
L, with respect the usual norm 

/ \ V2 

ll/(*,v)ll = V<7»7> = / Jf 2 (x,y)dxdy 

\ T / 

for / G £. 

Theorem 7.1. T/ie shortest element in L exists, and is given by 

n m mn 

where X(A) is the Lebesgue-measure of A. 

Proof. Since if is a closed linear subspace of C and L is just a shift of 
H, L has a shortest element fo(x,y) indeed. This fo(x,y) is uniquely 
determined by the following two properties: 

• fo(x,y) has column integrals Cf A (x) (x G [0, to]) and row inte- 
grals r fA (y) (y G [0, n\), 

• fo(x,y) is orthogonal to H, i.e. (fo(x,y), h(x,y)) = for every 
h(x,y) G H. 

We prove that the choice for fo(x,y) in the statement meets these 
requirements. To prove the first property, observe that 



r f A (y) d y = J J fA(x,y)dxdy — X(A). 

T 

Thus for any x G [0, to] we have 



J \ n to mn J 



, , X(A) X(A) 
(x) + ^-^ - n--^ = c fA (x). 
to mn J 



BOUNDS FOR DISCRETE TOMOGRAPHY SOLUTIONS 17 
Similarly, for any y G [0, n] 

" cjm + um _ m\ ^ m + rtM _ m m = r/M , 

n m mn J n mn 

o 

which proves the first property. 

In order to check the second property, take an arbitrary h(x, y) G H. 
For any x G [0, m] and y G [0, n] we have 

n m 

c h{x) = J h(x,y)dy = and r h (y) = J h(x,y)dx = 0. 



Then 



n m mn / 




h(x,y) 



o o 

m / n \ n 



\ n m mn J 



\0 / \0 / 

and the inner integrals are 0. This proves the second property, and the 
theorem follows. □ 
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